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OJJ' Abstract. We compute an approximate set of longitudinal quasiparticle modes for a 

3 ^ hydrogen condensate as produced in the MIT experiments. An expansion in quasiparticles 

■^^ , using a simple one-dimensional Bogoliubov picture shows however that at the high 

■ temperatures (f» 44/iK) and in the very shallow trap employed (uiz = 2tt X 10.2Hz) the 

O^ ' contribution to the density from the quasiparticles is about 20% of that from the condensate 

^\| ' mode, leading to an effective 32 (x, x) which varies between 1 and 3 depending on the position 

in the condensate. 
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1. Introduction 



cn 

^~^ i The analysis of the experimental data on a spin-polarized atomic hydrogen Bose-Einstein 

f^ I condensate as given in [^ g, ||, |||, relates the density of the condensate or of the non- 

^D . condensed gas to the frequency for the two-photon ls-2s transition by the cold collision 

line-shift formula 

A f 2ng2(x,x)nis 

m 



B 

' O ' In this formula ais-is and ais_2s are respectively the scattering lengths for collisions between 

C , the two atoms in the Is state or of one atom in each state, nis is the density of the Is atoms, 

and 



(7/^t(x)^t(xOV,(xO^(x)) 

52(x, X ) = — -— . (2) 

nis(x)nis(x') 

For a non-condensed thermal gas 52 (x, x) — 2, while for a pure condensate 52 (x, x) = 1. 
C^ ' The derivation of the formula dl]) is given in [b|1, though it is acknowledged there that 

this derivation is only valid for special cases; that is for a spatially homogeneous thermal 
non-condensed gas, and for a pure condensate. Intermediate situations, such as those which 
appear to pertain in practice have not been treated, and indeed are very difficult to analyze. 

In the experiments with the choice of 52 (x, x) = 1, frequency shifts measured in the 
Bose-Einstein condensate experiments would correspond to densities too high to be credible, 
and the data were analyzed using (72(x, x) = 2. This was justified by the manifest high 
temperature of the system, and the conjecture that condensation occuiTed into a multiplicity 
of states. 

This paper is the first of two on this subject. Here we shall concern ourselves specifically 
with the computation of 32 for the case of hydrogen condensate with a significant thermal 
component, in order to demonstrate what would be predicted by its use if the formula (p is 
assumed to be valid in that case. In our second paper [H] we shall show that the application of 
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^ to a condensate with a thermal component is not justified, and we shall give a full dynamic 
formulation of the actual processes involved. 

The evaluation of 52 is relevant to other phenomena, such as two-body loss processes, 
so its computation is still of value, even if it is not directly applicable to the analysis of 
the cold collision frequency shift in these experiments. Furthermore, in order to do the 
computations of our second paper, we shall need the eigenfunctions computed here. We 
will use a simple Bogoliubov model, though it will become apparent that there should still 
be significant corrections to this at the temperatures and densities found in practice. A more 
detailed study is left for a later paper 

One might expect that in general 32 (x, x) lies between 1 and 2, but in fact the only 
fundamental limit is that < .92 (x, x) < 00. However, in the case of Gaussian statistics, we 
get a stronger limit. Thus in the case of a fully thermalized gas, we can split the correlation 
function up using Gaussian factorization 

(^t(x)V.t(x')V(x')^(x)) = (^t(x)^(x))(^t(x')^(x')) 

+ (V/(x)^(x'))(Vt(x')VXx)) 
+ (Vt(x)Vt(x'))(V'(x')V(x)) 

(3) 
The third term in the case of a thermal gas is zero since there is no anomalous average 

(V^t(x)^t(x'))=0, (4) 

and setting x = x' gives 

, . _ (V't(x)^t(x)V,(x)V,(x)) 
^'^ ' ^ (^t(x)^(x))(^t(x)^(x)) • ^^^ 

Thus ^2 (x, x) = 2 is a consequence of the Gaussian nature of the fluctuations and of the 
absence of an anomalous average. The anomalous average does have an upper bound of 
approximately 

|(V'(x)^(x'))| < (^t(x)^(x)) (6) 

which gives the bound 

52(x,x)<3. (7) 

We will find in our calculation that in fact 52 (x, x) > 2 can occur in some regions of the 
Hydrogen vapor-condensate system for the experimentally realized parameters. 

2. Second-quantized Hamiltonian 

For the case under consideration the second-quantized Hamiltonian can be written in the form 



H= d 



j3. 



X 



{i^t(T + y(x))V' + ^V'WV'}- (8) 



i) The kinetic energy operator is T = —h V^ /2m 
ii) The interaction coefficient is given in terms of the scattering length by 

u — Anais^isfi /m (9) 

Experimentally 

ais_is = 0.0648nm (10) 

ais-2s = -1.4±0.3nm. (11) 
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3. Evaluation of 52 (x, x) for a warm condensate 

The hydrogen condensate at MIT is formed at 44/xK, and the radial and axial trap frequencies 
are 

w^ = 27r X 3.9kHz (12) 

^2 = 27r X 10.2Hz. (13) 

The ratio of these iOrjiOz — 382 gives the trap aspect ratio Izjlr, which at nearly 400:1 gives 
an almost one-dimensional condensate. 

As well as this, the ratios of the energies of the transverse and longitudinal quanta to the 
temperature are 

kTlhujr ~ 1/250 (14) 

fcr/?iw^ « 1/10^ (15) 

We can therefore expect that the average number of quasiparticles in the longitudinal mode 
is of the order of 10^ for the lower levels. This does not mean that the condensate is very 
different in size and shape, but its coherence may be affected by these large excitations. For 
example, the longitudinal Kohn mode corresponds to center of mass oscillations of the whole 
condensate, of mass lO^m, and using statistical mechanics, this means that 



V(^«3x 10"Vmw6x 10"^;^, (16) 

which is a completely unobservable deviation. 

3.7. Approximate evaluation of mode functions 

In order to get an approximate idea of the thermal fluctuations we will compute the 
wavefunctions of the lower lying excitations using a hydrodynamic approximation based 
on the work of Zaremba [|lO||. This method gives an approximate form for the fuU three 
dimensional wavefunction in the case that the radial modes are not excited; thus it is possible 
do a calculation by means of a one dimensional wave equation, and using the Thomas-Fermi 
approximation for the condensate wavefunction this give eigenfunctions analytically in terms 
of Jacobi polynomials. 

3.1.1. Hydrodynamic equations The system can be characterized by a density p(x, <) and a 
phase (/)(x, t), and the hydrodynamic equations are written in terms of the Unearized density 
fluctuation (5p(x, t) — p(x, t) — /Oo(x) as 

5p(x,i)= -Av.[po(x)V(/.(x,i)] (17) 

0(x,i) - --,5p(x,t). (18) 

Zaremba's analysis considers only the situation in which there is a cylindrically symmetric 
condensate with no z-dependence, and he shows that in this case the speed of sound along 
the z direction is half that expected in a homogeneous condensate with density equal to the 
peak density of the cylindrical condensate. We will generalize his result to the situation in 
which the trap is harmonic in all directions, but is very weak along the z-direction. Following 
Zaremba's lead, we look for the equation of motion for a perturbation with the property 

(5p(x, t) = bn{z) when po(x) 7^ 0, (19) 

= whenpo(x) ==0. (20) 
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Here it is assumed that the Thomas-Fermi form 



Po X - — -7> —. (21) 

is used, so that (|l9|) means that Sn{z) is only non-zero where there is a non-vanishing 
condensate density, i.e., when 



r < R{z) = J^ ^^. (22) 

When we substitute into the wave equation for 5p that arises from combining ( |l7| ) and([l8|), 
and then integrate over 27rr dr, the resulting equation can be written in terms of the scaled 
length 



z^hJ^^ = Xh (23) 



as 

:)2 



^t Ah \ 



-f{i^-h')^^~^h^\- (24) 



e ] 
eigenvalue spectrum 



The operator on the right hand side is that of the Jacobi polynomials P„ ' (ft,), leading to the 



LOn = -^yjn{7l + i). (25) 

The more exact three dimensional analysis of Fliesser et al. ^ gives, in the case of an 
extended cigar-shaped condensate, the same formula. 

Note that these Jacobi polynomials are orthogonal in the sense that 

[2n + 3)(ri + 2) 

(26) 

so that the appropriate weight function is (1 — h?), which is proportional to the stationary 
condensate density. 

3.1.2. Three dimensional interpretation of the wavefunctions As in Zaremba's analysis, 
these eigenfunctions are constant in the radial direction to the edge of the Thomas-Fermi 
condensate (Pl|), where they abruptly drop to zero. This means that the excitations are all 
essentially sound waves travelling along the length of the condensate. 

3.1.3. Quasiparticle wavefunctions We can thus write 

5n{z, t)=Yl ^« cosKt)P,(i^i) (h) (27) 

n 

4>{z,t) = -Y^^sin{uj^t)Pi^^^\h). (28) 

n 

The resulting condensate wavefunction ^ (normalized to the total number of particles, not to 
1) is given by 



^(x, t) = v/po(x)+^p(x,i)exp(i0(z, t)) (29) 



Analysis of g2 in hydrogen. . . 5 

<5p(x, t) 



Vpo(x) + -^^M= + Vpo(x) 4>{z, t) (30) 

Vpo(x) 

'1 upo(x) 






(31) 



2^po(x) V2 ?i^n 

This means that the equivalent u„ and u„ functions which would turn up in a quasiparticle 
expansion of the field operator, 

V'(x,i) = ^1 v/M^+^{«„e-^'^"*M„(x) +4e'""S„(x)} 1 (32) 

{^/M^ + Vi(x)} (33) 



...M^^^s^fi+^^je-'c) 



Oq 



/n. 
are given by 



2VPo(x) 



2VPo(x) 
where ^„ are normalization constants to be determined 



(35) 



3.1.4. Amplitudes orthogonal to the condensate The correct eigenfunctions whose quantized 
amplitudes represent quasiparticles are orthogonal to the condensate wavefunction. To check 
this, we evaluate 

d^x VPo(x)wn(x) = ^ / dh j 27rrdr v/po(x) 






Q,H^)p,u,„, 



/.-f'-K^^)^-'"'"* 



2 . . _ 

(37) 

We now use the Thomas-Fermi form ( pi| ) for the stationary condensate density so that 



/ 2TTr dr ^ 'KR{zf 

Jo 



(38) 



27rrdrMpo(x) = -^(1 - /(,2)2. (39) 







so that, using the orthogonality property (|6|), equation ( p7| ) becomes 
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The nonorthogonal part — which turns up only in the quadrupole mode — represents an 
unobservable time dependent phase of the condensate, and is correctly treated by simply 
subtracting the component parallel to the condensate [Bl 0], equivalent to making the 
replacement in the second part of the large bracket in ( |34| ) and ( p^ P„ ' {h) ^ Pn (h), 
where 

Pi'''Hh) = Pi'''Hh) n^2, (41) 

P^'''Hh)^^^{7h'~l), (42) 

where the norm of P2 ' has been fixed to be the same as that of Pj ^ namely 6/7 

It is als( 
other, so that 

d^x[u„(x)w„j(x) - w„(x)w„j(x)] 



It is also straightforward to check that all the P„ ' for n > 1 are orthogonal to each 



/' 



jniiuAl 

= A ^- On,niR-n (43) 



where 




^"" "^ (27i + 3)(n + 2)+^"'' 



/24 6 
35 ^7 



(44) 



Since the normaUzation is given by 

/ d^x[u„(x)w„(x) - ■y„(x)w„(x)] = 1 (45) 

we must choose 



3.7.5. r/ie number of particles per quasiparticle Neglecting the zero point contribution, 
negligible in comparison to the thermal contribution in this case, the number of particles per 
quasiparticle for the mode n is given by 

d3x[u„(x)2+i;„(x)2] (47) 



= Ai d 



j3. 



x 



P^^^^\hf , u'po{^)Pk^'^\hf 



8/30 (x) 2^2 w2 j ■ 

(48) 
The first term involves a divergent integral, arising from the form (Eni for po (x) . This problem 



is a reflection of the failure of the linearized expansion (pOh This can be fixed by noting that the 
true density does not go to zero at /i = 1, but reaches a value dependent on the healing length. 
The result is only logarithmically divergent, and any estimate of it is very much smaller than 
the second term, which does not diverge. We will therefore neglect this first term. 
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Figure 1. Comparison of atom numbers in the lowest quasiparticle modes for a sodium and 
a hydrogen condensate, a) Solid line: occupation of each mode; b) Dotted line: cumulative 
number up to the given level. 



Using the expression ( pi| ) for /5o(x), we need the results 

dh{l-h''f[P^^^^\h)Y=Sr, 

8(n+l) 8(n + l) ( {n+\){n + 'i) r 





(49) 


(n^2) 


(50) 


{n ^ 2) 


(51) 



(2n + 3)(n + 2) (2n + 3)2 \ (n + 2)(2n + 5) [In + 1) J ' 



21' 
and using the expHcit form for A„ (p6[), we find 

2huJn'Ti-n 

Here p — fi/u is the peak condensate density, which in the experiment of Fried et al. 
reported to be 4.8 x lO^^m"'^. 



(52) 



3.2. Numbers and densities for the experimental situation 

Using the resuh ( P5| ) and the experimental value for lOz, we have w„ = 27r x 10.2Hz x 
y^n{n + 3)/4, leading to Qn ~ 37765„/-\/n(ri + 3)7^„. This would lead to 



Nn 



kT _ 6.9 X IQ^Sn 
hujn n(n + 3)TZn 



(53) 
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Figure 2. Plots showing: a) The individual confributions to the total density from the 
condensate and the first 10 eigenfunctions; b) From bottom to top, the condensate density, the 
total density, and the total density multipUed by 92 (x, x) (Normalization for both chosen so 
that the peak density arising from the condensate wavefunction is 1); c) The 512 (x, x) arising 
from all of the components. 

particles in each mode. In Fig.|l]we illustrate the contributions to the particle numbers from 
the condensate and the first two condensate modes for the hydrogen condensate and for a 
typical MIT sodium condensate [^, whose geometry is a much less extreme cigar shape. The 
net contribution from the quasiparticle modes is about 2% of the number in the condensate 
mode for sodium, but for Hydrogen approaches 20% of the number in the condensate mode. 
Notice that convergence appears after about 10 modes. 



3.2.1. Particle density arising from quasiparticles The contribution to n{h) arising from the 
quasiparticles is simply of the form of the integrand in (EtI) for each mode, with the resulting 
overall density taking the form 



"qW = Yl 



3(14) 



AlkT PrV'^'ik)^ , u^poi^)Pk'''>{h) 



(1,1), 



hujn [ 8/90 (x) 2?i^w2 

We wiU need the anomalous average as well, that is the term 

and this is quite readily shown to be 



T-Q 



w = E 



AlkT i Pi'^'\h)^ u^p,{^)Pi'^'\h)^ 



huJn 



8po(x) 



27i w. 



2,„2 



(54) 



(55) 



(56) 
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Figure 3. Plots vs. r and z showing: a) The condensate density (in the Thomas-Fermi 
approximation); b) The density arising from the condensate plus first 10 quasiparticle states 
(thermally occupied); c) The effective density, obtained from b) by multiplying by 32 (x, x) 



In the degree of approximation we are using, in which the first terms in the bracketed parts of 
the right hand sides of (gj) and ( p6| ) are neglected, it follows that mQ{x) = —iiqIx.) — the 
anomalous term is thus very important. 

3.2.2. Computation ofg2{x, x) The quasiparticle contribution to the field operator 

^(x) = Y. {«ne^^""*w„(x) + 4e*""*i;„(x)} (57) 

n 

is composed of operators a„, aj^, for which the density operator is quantum Gaussian with 
zero mean, so the Gaussian factorization of the the four point correlation function is 

Jv3t(x)^t(x)v:(x)V3(x)) = 2nQ(x)2 + |mQ(x)2| (58) 

« 371q(x)2. (59) 

Thus, for the full field operator 

(V;t(x)Vt(x)V;(x)^(x)) = po(x)2 + 4po(x)nQ(x) 
+ 2po(x)toq(x) + 2nQ(x)^ + |mQ(x)2| 
w po(x)2 + 2po(x)nQ(x) + 3nQ(x)2 
Noting now that the particle density is 
n(x) =: po(x) +nQ{x.) 
= Po(x)[l + <z(x)] 
which defines g(x), we see that 

Notice that the possible values of 172 (x, x) are between 1 and 3, and that 52 (x, x) does depend 
on position. 

Using these expressions we can now compute both n{x.) and (72(x, x), and these are 
illustrated in Fig.0. 



(60) 
(61) 

(62) 
(63) 



(64) 
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3.3. Interpretation 

Fig.0 consists of three plots, to illustrate the two major effects arising from this calculation. 
Comparing a) with b), it can be seen how the condensate density is modified both in shape and 
size by the excitation of the quasiparticles. Comparison of b) with c) shows how inclusion of 
the factor of 52 (x, x) would modify the signal from the experiment not by simply multiplying 
by a factor of 2 as has been assumed in the experiment. 

4. Conclusion 

The results of this calculation are only of indicative interest for the following reasons: 

a) The occupations of the quasiparticle levels are large enough that the simple Bogoliubov 
formulation may need modification. 

b) There is as yet no justification for simply including the factor of 52(x,x) as the 
appropriate correction for correlations. 

c) The 92 (x, x) computed is the result of an effective interpolation (given by (^J)) of 
52 (x, x) ~ 3 for the quasiparticles and (72 (x, x) ^ 1 for the condensate. It can be shown 
in fact that the experimental frequency shift which results in the hydrogen condensate 
experiments is different for condensate and quasiparticles, so that there is no simple 
interpretation in terms of a single frequency shift determined by 32 (x, x). 

These matters are be attended to in our second paper [0], but it is already clear that the correct 
interpretation of the cold-collision shift in these experiments is a matter of great significance. 
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